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Abstract 

We study the maximal displacement of branching random walks in a class of time inhomoge- 
neous environments. Specifically, binary branching random walks with Gaussian increments will 
be considered, where the variances of the increments change over time macroscopically. We find 
the asymptotics of the maximum up to an Op(l) (stochastically bounded) error, and focus on the 
following phenomena: the profile of the variance matters, both to the leading (velocity) term and 
to the logarithmic correction term, and the latter exhibits a phase transition. 

1 Introduction 

Branching random walks and their maxima have been studied mostly in space-time homoge- 
neous environments (deterministic or random). For work on the deterministic homogeneous 
case of relevance to our study we refer to |1] and the recent [1] and [2]. For the random 
environment case, a sample of relevant papers is [121 [El El HI 1211 [221 [22]. As is well 
documented in these references, under reasonable hypotheses, in the homogeneous case the 
maximum grows linearly, with a logarithmic correction, and is tight around its median. 

Branching random walks are also studied under some space inhomogeneous environments. 
A sample of those papers are [71 [9l [IH US HH [I9] . 

Recently, Bramson and Zeitouni [5] and Fang [10] showed that the maxima of branching 
random walks, recentered around their median, are still tight in time inhomogeneous envi- 
ronments satisfying certain uniform regularity assumptions, in particular, the laws of the 
increments can vary with respect to time and the walks may have some local dependence. 
A natural question is to ask, in that situation, what is the asymptotic behavior of the max- 
ima. Similar questions were discussed in the context of branching Brownian motion using 
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PDE techniques, see e.g. Nolen and Ryzhik [21], using the fact that the distributions of the 
maxima satisfy the KPP equation whose solution exhibits a travehng wave phenomenon. 

In all these models, while the linear traveling speed of the maxima is a relatively easy 
consequence of the large deviation principle, the evaluation of the second order correction 
term, like the ones in Bramson [1] and Addario-Berry and Reed [1], is more involved and 
requires a detailed analysis of the walks; to our knowledge, it has so far only been performed 
in the time homogeneous case. 

Our goal is to start exploring the time inhomogeneous setup. As we will detail below, 
the situation, even in the simplest setting, is complex and, for example, the order in which 
inhomogeneity presents itself matters, both in the leading term and in the correction term. 

In this paper, in order to best describe the phenomenon discussed above, we focus on the 
simplest case of binary branching random walks where the diffusivity of the particles takes 
two distinct values as a function of time. 

We now describe the setup in detail. For a > 0, let A^(0,cr^) denote the normal distri- 
butions with mean zero and variance . Let n be an integer, and let (jj,cr| > be given. 
We start the system with one particle at location at time 0. Suppose that f is a particle 
at location at time k. Then v dies at time /c + 1 and gives birth to two particles v\ 
and f 2, and each of the two offspring ({f i = 1, 2}) moves independently to a new location 
Sm with the increment Sm — independent of Sy and distributed as N{0,al) if k < n/2 
and as A^(0, cxg) if ^/2 < k < n. Let D„ denote the collection of all particles at time n. 
For a particle v G D„ and i < n, we let denote the ith level ancestor of v, that is the 
unique element of Dj on the geodesic connecting v and the root. We study the maximal 
displacement Mn = maXt,giD„ Sy at time n, for n large. 

It will be clear that the analysis extends to a wide class of inhomogeneities with finitely 
many values and 'macroscopic' change (similar to the description in the previous paragraph), 
and to the Galton- Watson setup. A universal result that will allow for continuous change of 
the variances is more complicated, is expected to present different correction terms, and is 
the subject of further study. 

In order to describe the results in a concise way, we recall the notation Op{l) for stochas- 
tically boundedness. That is, if a sequence of random variables Rn satisfies Rn = Op(l), 
then, for any e > 0, there exists an M such that P(|-R„| > M) < e for all n. 

An interesting feature of M„ is that the asymptotic behavior depends on the order relation 
between af and ctj. That is, while 

M„ = (^21^ a^s) n - /3^|= ^ogn + Op(l) (1) 

is true for some choice of Ooff and /3, cigg and (3 take different expressions for different ordering 
of 0"! and (72 • Note that Q is equivalent to say that the sequence {M„ — Med{Mn)}n is tight 
and 

Med{M^) = (v/21og2aeff)n-/3 ^^^^ logn + 0(1), 

^ Since one can understand a branching random walk as a 'competition' between branching and random 
walk, one may get similar results by fixing the variance and changing the branching rate with respect to 
time. 
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where Med{X) = sup{x : P{X < x) < ^} is the median of the random variable X. In the 
following, we will use superscripts to distinguish different cases, see (g, (g and Q below. 

A special and well-known case is when cxi = (72 = ex, i.e., all the increments are i.i.d.. In 
that case, the maximal displacement is described as follows: 

M= = (^21^ a) n - ^-^^=logn + Op(l); (2) 

the proof can be found in |T], and its analog for branching Brownian motion can be found 
in |1] using probabilistic techniques and [20] using PDE techniques. This homogeneous case 
corresponds to ([T]) with cTgg = a and /9 = |. In this paper, we deal with the extension to 
the inhomogeneous case. The main results are the following two theorems. 

Theorem 1. When a\ < (j| (increasing variances), the maximal displacement is 



Ml = {J{al + al) log 2^ n - logn + Op(l), (3) 



which is of the form ([T| with a^g = y "^2"^ ^'^^ P ~ \- 

Theorem 2. When af > cr| (decreasing variances), the maximal displacement is 

Mi = n - ^-^== log n + Op(l), (4) 

which is of the form ([T| with = ^i±^ and /3 = 3. 

For comparison purpose, it is useful to introduce the model of 2" independent (inhomo- 
geneous) random walks with centered independent Gaussian variables, with variance profile 
as above. Denote by M™*^ the maximal displacement at time n in this model. Because of 
the complete independence, it can be easily shown that 



i 



Ml^'' = ( Jicrj + al) log2 ) n - l^gn + Op(l) (5) 



for all choices of af and cxj. Thus, in this case, cXgg = 

Vi^l + ^i)/2 and /3 = 1/2. Thus, 
the difference between and M™'^ when = (j| lies in the logarithmic correction. As 
commented (for branching Brownian motion) in [3], the different correction is due to the 
intrinsic dependence between particles coming from the branching structure in branching 
random walks. 

Another related quantity is the sub-maximum obtained by a greedy algorithm, which only 
considers the maximum over all decendents of the maximal particle at time n/2. Applying 
(|2|, we find that the output of such algorithm is 

n 3 (Ti T f r~-\ n '6 02 , 



^^^^^^^^2 - 27^^°^ 2 J + ^v/21og2a2- - ^ 7^ 2 J + ^^^^^ 
y2bp(ai + (T2) 3(01 + 02), ^ „ 
2 " - ^ 
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Comparing ([6]) with the theorems, we see that this algorithm yields the maximum up to an 
Op{l) error in the case of decreasing variances (compare with Q) but not in the case of 
increasing variances (compare with ([3])) or of homogeneous increments (compare with ([2])). 
A few comparisons are now in order. 

1. When the variances are increasing, is asymptotically (up to Op (1) error) the same as 
M^'^, which is exactly the same as the maximum of independent homogeneous random 
walks with effective variance 

2. When the variances are decreasing, M;^ shares the same asymptotic behavior with the 
sub-maximum In this case, a greedy strategy yields the approximate maximum. 

3. With the same set of diffusivity constants {cr^,cr|} but different order, is greater 
than Mj^. 

4. While the leading order terms in ([2]), (|3]) and ^ are continuous in ai and a2 (they 
coincide upon setting ai = (72), the logarithmic corrections exhibit a phase transition 
phenomenon (they are not the same when we let (Xi = (J2). 

We will prove Theorem [T] in Section [2] and Theorem |2] in Section |3j Before proving the 
theorems, we state a tightness result. 

Lemma 1. The sequences {Ml - Mec?(M^)}„ and {M^ - Med(M^)}„ are tight. 

This lemma follows from Bramson and Zeitouni [5] or Fang pl)j. One can write down 
a similar recursion for the distribution of M„ to the one in those two papers, except for 
different subscripts and superscripts. Since the argument there depends only on one step of 
the recursion, it applies here directly without any change and leads to the tightness result 
in the lemma. 

A note on notation: throughout, we use C to denote a generic positive constant, possibly 
depending on ai and (72, that may change from line to line. 



2 Increasing Variances: af < (j\ 

In this section, we prove Theorem [T| We begin in Subsection 2J_ with a result on the 



fluctuation of an inhomogeneous random walk. In the short Subsection 2.2 we provide large- 
deviations based heuristics for our results. While it is not used in the actual proof, it explains 
the leading term of the maximal displacement and gives hints about the derivation of the 



logarithmic correction term. The actual proof of Theorem [T] is provided in subsection 2^ 



2.1 Fluctuation of an Inhomogeneous Random Walk 

Let 

n/2 n 
2=1 i=n/2+l 
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be an inhomogeneous random walk, where ~ A^(0, cr^), ~ A^(0, erf), and and are 
independent. Define 



ggf + a|(fc-f) n 



n 
2' 



< k < n, 



(8) 



and 



k.n 



(9) 



A; < n/2, 

Cf{n-kf/^, n/2<k<n. 

As the foUowing lemma says, conditioned on {Sn = x}, the path of the walk Sn follows 
Sk,n{x) with fluctuation less than or equal to fk,n at level k <n. 

Lemma 2. There exists a constant C > (independent of n) such that 

,nt Sk.n {Sn) + fk,n] for all 0<k< n\Sn) > C, 

where Sn{k) is the sum of the first k summands of Sn, i-e., 

( k 

A;<n/2, 



Sn{k) 



k=l 

n/2 



fe=n/2+l 



n. 



k=l 



Proof. Let Sk^n = Sn{k) — Sfc,n('S'n)- Then, similar to Brownian bridge, one can check that 
Sk,n are independent of Sn- To see this, first note that the covariance between Sk^n and Sn is 

,m Sn) — ESk ,nSn — ESk^nESn — ESk ,nSni 



since ESn — and ESk,n — 0. 
For k < n/2. 



Sk,; 



1 - 



alk 



n/2 



alk 



2/ i=l 



al + al) 



2 i=k+l 



(al + al) 



2 i=n/2+l 



Expand Sk,nSn, take expectation, and then all terms vanish except for those containing Xf 
and Y-^. Taking into account that EXf — al and EY^ — a 2, one has 



Cov{Sk,n, Sn) — ESk,nSn 



= 1 



alk 



2/ i=i 



kal- 



alk 



2U '^/^ 

2^Y.EX^ 

i=k+l 



[ai + a^j^ 



alk 



[a^ + a2j2 



=n/2+l 



(10) 
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For n/2 < k < n, one can calculate Cov{Sk,n, Sn) = similarly as follows. First, 
~ _ aj{n - k) ^ y <yl{n-k) A ^ _ A _ <yl{n - k) \ >A 

Then, expanding Sk^nSn and taking expectation, one has 

Cov{Sk,ni Sn) = ESk,nSn 

K^l +^2)2 i=i 1^1 +^2^2 ,=n/2+l ^ 1^1 + ^2J 2 / 
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= 

Therefore, Sk,n are independent of Sn since they are Gaussian. Using this independence, 

P{Sn{k) e [Sk,n{Sn) " fk ,ni Sk.n {Sn) + fk,n] for all < k < n\Sn) 
= P{Sk,n e [-fk,n, fk,n] foT all < A; < n\Sn) 
= PiSk,n e [-fk,n, fk,n] for all < A; < n). 



By calculation similar to (|10|), Sk^n is a Gaussian sequence with mean zero and variance 

2 

quantity is 



2.?.) ^ ^ ^ _ (K^+. )„ 2.i(n-.)) < A: < n. The above 



1 - Pi\Sk,n\ > fk,n, for some < k < n) > I - "^P^Sk^nl > fk,n)- 

k=l 

Using a standard Gaussian estimate, e.g. [Bl Theorem 1.4], the above quantity is at least, 

1 - y ^e^%- > 1 - 2 y ^e-/-'^^^^ -.= 

where cq, ci are constants depending on ai and (T2, and C > can be realized by choosing 
the constant Cf large. This proves the lemma. □ 



2.2 Sample Path Large Deviation Heuristics 

We explain (without giving a proof) what we expect for the order n term of Mnt, by giving 
a large deviation argument. The exact proof will be postponed to the next subsection. 
Consider the same Sn as defined in ^ and a function </)(t) defined on [0, 1] with 0(0) = 0. 
A sample path large deviation result, see |6J, Theorem 5.1.2], tells us that the probability for 
S^rn} to be roughly (f){r)n for < r < s < 1 is roughly e~'^^'''^'^\ where 

J,(0)= rA*Mr))dr, (11) 
^0 
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dr *" 



and A*(x) 



2 2 

for < r < 1/2, and f-2, for 1/2 < r < 1. A first moment 



argument would yield a necessary condition for a walk that roughly follows the path 0(r)n 
to exist among the branching random walks, 



< slog2, for all < s < 1. 



:i2i 



This is equivalent to 



20-2 



-dr < slog 2, 



2a( 



-dr 



2al 



-dr < slog 2, 



< s < 



< s < 1. 



(13) 



Otherwise, if (12) is violated for some sq, i.e. 
(j){r)n to 



(0) > So log 2, there will be no path following 



(so)?T., since the expected number of such paths is 2*"e 
which decreases exponentially. 



-nls{<l>) 



- (I sW-s log 2)n 



Our goal is then to maximize 0(1) under the constraints (13). By Jensen's inequality 



and convexity, one can prove that it is equivalent to maximizing 0(1) subject to 
<^'(l/2) 0^(1/2) , (0(l)-0(l/2))2 



at 



1, 0'(l/2) 



< log 2. 



(14) 



Note that the above argument does not necessarily require < erg- 

Under the assumption that af < cr|, we can solve the optimization problem with the 
optimal curve 



0(s) 



2a?v45g2 1 2cT|v^bp 



1 



a. 



2 



+ 



(To 



< s < 



< s < 1. 



(15) 



If we plot this optimal curve and the suboptimal curve leading to ^ as in Figure [T| 
it is easy to see that the ancestor at time n/2 of the actual maximum at time n is not a 



maximum at time n/2, since 



2(TjVlog2 



< ^2a\ log 2. A further rigorous calculation as in 



the next subsection shows that, along the optimal curve (15), the branching random walks 



have an exponential decay of correlation. Thus a fluctuation between n}!"^ and n that is 
larger than the typical fluctuation of a random walk is admissible. This is consistent with 
the naive observation from Figure fTl This kind of behavior also occurs in the independent 

— ' A 

random walks model, explaining why M\ and have the same asymptotical expansion 

up to an 0(1) error, see (§ and (§. 

2.3 Proof of Theorem [I] 



With Lemma |2] and the observation from Section 2f2_. we can now provide a proof of Theorem 
[T| applying the first and second moments method to the appropriate sets. In the proof, we 
use Sn to denote the walk defined by ([T]) and S'fc to denote the sum of the first k summand 
in Sn- 
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Figure 1: Dashed: maximum at time n of BRW starting from maximum at time n/2. 
Solid: maximum at time n of BRW starting from time 0. 



Proof of Theorem\^ Upper hound. Let = + cr|) log2n— ^^J^^ logn. Let N\^n = 

^yeD„ '^{Sv>a„+y} be the number of particles at time n whose displacements are greater than 
a„ + y. Then 

EiVi,„ = 2"P(5, >an + y)< cae"^^^ 

where C2 and C3 are constants independent of n and the last inequality is due to the fact that 
Sn ~ ^(0, ^^^^^n). So we have, by the Chebyshev's inequality, 

P{Ml >an + y) = P{Ni > 1) < EiVi,„ < cae-^^^. (16) 

Therefore, this probability can be made as small as we wish by choosing a large y. 

Lower bound. Consider the walks which are at Sn & In = [o-n, cin + 1] at time n and follow 
Sfc,n(sn), defined by ([s]), at intermediate times with fluctuation bounded by fk^m defined by 
(§. Let h, (x) = [sk^n{x) — fk,n, Sk,n{x) + fk,n] bc the 'admissible' interval at time k given 
Sn = X, and let 

N2,n = ^ -'-{s„e/„,5^fce/fc,„(s„) for all o<fc<n} 
be the number of such walks. By Lemma [2| 

EN2,n = 2''P{SneIn,Sn{k)eh,n{Sn) for al\0<k<n) 

= 2"E(1|5„6/„}P(5„(A;) G IkASn) for all < A; < n\Sn)) 

> 2"CP(5„ G 4) > C4. (17) 

Next, we bound the second moment EN^^. By considering the location of any pair 
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Vi,V2 e D„ of particles at time n and at their common ancestor Vi Av2, we have 
E^2,n = E XI ^{Sy^eI„,s, .jeij,n{s,^.)j) ior all o<j<n,i=i,2} 



{Svi^in, s^^,~^jeij,„{s^^,s^j) for all o<j<n,i— 1,2} 



k=0 fi,f2SD„ 



< XI 1] ^^'^'^1 ^ '^('^i)^ ^ ^j,niS(v,y) for all < J < n) 



fc=0 Vl,V2eT)n 



'P{S^,,^ 'S'ujAy2 ^ [-^ Sk,n{x) fk,m^ Sk,n{x) +/fc,n])2^ ^ -^n); 

where we use the independence between 8^2 — 5*^,^^1)2 ^{viY the last inequality. And 
the last expression (double sum) in the above display is 



X2'"-'P(^n e In:Sn{j) G IjA^n) for all < J < n) 

k=0 

'P{Sn ~ Sn{k) E [x — Sk,n{x) — fk,m X — + fk,n]j X G /„) 

n 

< EN2,n X '^''~^P{Sn - Sn{k) G [x - Sk,n{x) - fk,n, X - Sk,n{x) + fk,n],X G Q ■ 



k=0 



The above probabilities can be estimated separately when k < n/2 and n/2 < k < n. For 
k < n/2, Sn - Sn(k) ~ A^(0, f (a? + a^) - kaf). Thus, 



P{Sn Sn{k) G [a; SA;,n(^) fk,mX Sk^n{x) + /fc^n]?^; G /„ 



< 2/; 



k.n 



^TTiiaf + ai)n - 2kal) 



exp 



/ 



V 



((1 (o-2+a2)n)'^n /m) 



< 2 '^l+''2 

For n/2<k< n, Sn - Snik) ~ iV(0, (n - A;)(7|). Thus, 

P{Sn ~ Sn{k) G [x — Sk,n{x) — fk,n, X — Sk,n{x) + fk,n]i X G In) 

( (2al(n-k)^ f_ \^\ 



^27r(n — fc)cr; 



: exp 



v 



2a|(«-fc) _ f V 

2{n - k)al 



J 



^ cT (n-fc)+o(n-fc) 



Therefore, 



"M / 2 2 2 

EiVj„<£7V,„|$:2^F3'«'"+ 5: 2 



fc=0 



_ 1+^2 

A;=n/2+l 



i '^(n-k)+o{n-k) 



< C5EN2,n, (18) 



k+o{k) 



where C5 = 2^^g2'^i+'^2 . By the Cauchy-Schwartz inequahty, 

P{Ml > a„) > P(iV,,„ > 0) > > C4/C5 > 0. (19) 



The upper bound (16) and lower bound (19) imply that there exists a large enough 
constant yo such that 

P{Ml G [an,an + yo]) > ^ > 0. 

ZC5 

Lemma [1] tells us that the sequence {M^ — Med(M^)}„ is tight, so = a„ + 0(1) a.s.. 
That completes the proof. □ 

3 Decreasing Variances: af > (j\ 

We will again separate the proof of Theorem |2] into two parts, the lower bound and the 
upper bound. Fortunately, we can apply ([2]) directly to get a lower bound so that we can 
avoid repeating the second moment argument. However, we do need to reproduce (the first 
moment argument) part of the proof of ([2]) in order to get an upper bound. 

3.1 An Estimate for Brownian Bridge 

Toward this end, we need the following analog of Bramson [U Proposition 1']. The original 
proof in Bramson's used the Gaussian density and reflection principle of continuous time 
Brownian motion, which also hold for the discrete time version. The proof extends without 
much effort to yield the following estimate for the Brownian bridge — \Bni where is 
a random walk with standard normal increments. 

Lemma 3. Let 

{0 if s = 0,n, 

100 log A; z/A; = l,...,n/2, 

1001og(n-A;) ifk = n/2,...,n-l. 

Then, there exists a constant C such that, for all y > 0, 

k C(l + 

P{Bn - -Bn < L{k) +y forO<k<n)< ^ ^ ^> 



n n 

The coefficient 100 before log is chosen large enough to be suitable for later use, and is 
not crucial in Lemma |3l 



3.2 Proof of Theorem H 



Before proving the theorem, we discuss the equivalent optimization problems ( 13 ) and ( 14 ) 
under our current setting cr^ > cr|. It can be solved by employing the optimal curve 

V21og2ais, < s < J, 

(20) 

V2 log 2^1^ + V21og2a2(s - ^), ^ < ^ < 1- 
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If we plot the curve 0(s) and the suboptimal curve leading to ^ as in Figure |2| these two 
curves coincide with each other up to order n. Figure [2] seems to indicate that the maximum 
at time n for the branching random walk starting from time comes from the maximum at 
time n/2. As will be shown rigorously, if a particle at time n/2 is left significantly behind the 
maximum, its descendents will not be able to catch up by time n. The difference between 
Figure [l] and Figure [2] explains the difference in the logarithmic correction between and 

n 




Figure 2: Dash: both the optimal path to the maximum at time n and the path leading 
to the maximum of BRW starting from the maximum at time n/2. Solid: the path to the 
maximal (rightmost) descendent of a particle at time n/2 that is significantly behind the 
maximum then. 



Proof of Theorem\^ Lower Bound. For each i = 1,2, the formula ^ implies that there 
exist Hi (possibly negative) such that, for branching random walk at time n/2 with variance 



n 3(T,- , n 



P ( M„,, > ^2 log 2a.- - log 2 + !/. j > 2- 



By considering a branching random walk starting from a particle at time n/2, whose location 
is greater than a/2 Ic 
and 2, we know that 



is greater than a/2 log 2cri| — 2^2'iog2 ^ + ^"^^ applying the above display with i = 1 



( I J2 log2(ai + 02) 3(ai + (12) , ^ \ 1 

P (Mi > ' '' n - log - + j > (21) 

Upper Bound. We will use a first moment argument to prove that there exists a constant 
y (large enough) such that 



/ ^21og2(..+.,) _ 3(.i+...) n^\^^ 

V 2 2^21^ ^2 10 ^ ^ 
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\/2T5f2(ai+f72) 



Similarly to the last argument in the proof of Theorem [T| the upper bound (22) and the 
lower bound (21), together with the tightness result from Lemma[l| prove Theorem[2] So it 
remains to show (22). 

Toward this end, we define a polygonal line (piecewise linear curve) leading to 
^^=i log f as follows: for 1 < < n/2, 

M{k) 

and for n/2 + 1 < k < n, 



n— 



k I— n 3(Ti 

— (v/21og2a,--^-^log-J; 



n 



3(79 



M(k) = M(„/2) + ^( v/2tog2..- - 



Note that - logn < log k for k < n. Also define 



k 0, 2 5 ^5 



log 



fik) 



y 








y + 


5(71 


log k 




2V2 log 2 




y + 




log(f 


-k) 


2^2 log 2 


y + 


5(72 


log (A; 




2V2 log 2 


-f) 


y + 


5(72 


log(n 


-k) 


2V2 log 2 



l<k< n/4, 
A;) f < A; < f - 



3n 



3n 



< A; < n - 1. 



We will use f{k) to denote the allowed offset (deviation) from M{k) in the following argu- 
ment. 

The probability on the left side of (22) is equal to 

P{3v G D„ such that Sy > M{n) + y). 
For each v G 1D„, we define = inf{A; : S^k > M{k) + /(A;)}; then (22) is implied by 

n 

P{^v G D„ such that S'„ > M{n) +y,Ty = k) < 1/10. (23) 

k=l 

We will split the sum into four regimes: [l,r;,/4], [n/A,n/2], [n/2,3n/A] and [3n/A,n\, corre- 
sponding to the four parts of the definition of f{k). The sum over each regime, corresponding 
to the events in the four pictures in Figure [3| can be made small. The first two are the dis- 
crete analog of the upper bound argument in Bramson [4j . We will present a complete proof 
for the first two cases, since the argument is not too long and the argument (not only the 
result) is used in the latter two cases. 

(i). When l<A;<n/4, we have, by the Chebyshev's inequality, 

P{3v G D„ such that > M{n) +y,Ty = k) 
< Pi3v G Bk, such that Sy > M{k) + f{k)) < E ks.>M{k)+m}- 
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(a) 1 to n/4 (b) n/4 to n/2 




(c) n/2 to 3n/4 (d) 3n/4 to n 



Figure 3: Four small probability events. Dash line: M{k). Solid curve: M{k) + f{k). 
Polygonal line: a random walk. 



The above expectation is less than or equal to 



C2 



k (M(k)+S{k)Y 



< 



CT 



exp 



V21og 2(71 A; 



s/2 log 2 



log /c + y 



v 



2kal 



Summing these upper bounds over k e [l,n/4], we obtain that 



n/4 



P{3v e D„ such that > M{n) + y,Ty ^ k) <Ce ^ ^ A;" 



-3/2 



(24) 



(25) 



fc=i 



ifc=i 



The right side of the above inequality can be made as small as we wish, say at most j^, by 
choosing y large enough. 

(ii). When n/A < k < n/2, we again have, by Chebyshev's inequality, 

P{3v G Dn such that > M{n) +y,T^ = k) 
< P{3v e Dfc, such that Sy > M{k) + f{k), and S^i < M{i) + f{i) ior 1 < i < k) 



< 



{Sv>M{k)+f{k), and s„i<M(i)+/(i) for i<i<fe}" 



Letting Sk be a copy of the random walks before time n/2, then the above expectation is 
equal to 

2''P{Sk > M{k) + f{k), and < M{i) + f{i) for 1 < i < A;) 
< 2^P{Sk > M{k) + f{k), and -{S^ - '-Sk) < -{f{i) - Ifik)) for 1 < i < k). 

(26) 
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{Si — j:Sk) is a. discrete Brownian bridge and is independent of Sk- Because of this inde- 



pendence, the above quantity is less than or equal to 

2^PiSk > M{k) + f{k)) ■ P{-{Si - jSk) < -(fit) - jfm for 1 < 2 < k). 



The first probability can be estimated similarly to ([24j), 

P{Sk > M{k) + f{k)) 



C 



V 



2kal 



< C2-'kCl-k)-'l\-'W\ 



(27) 



To estimate the second probability, we first estimate ^{f{i) — ^f{k)). It is less than or 

equal to ^f{i) = ^ + 2 V2 2 for i ^ < n/A, and, for k/2 < i < A;, it is less than 
or equal to 

log(n/2 - i) - log(n/2 -k) + log(n/2 - A;) ) + — (1 - ' 



< 



2y2bp 
5 

2v/2bP 



log(fc + ^^-r-^ log ) + — < 100 log(A; - i) + 



k 



0"! 



Therefore, applying Lemma |3j we have 



P (-{S. - Is,) < -(/(z) - lf{k)) for 1 < z < A: 
\ai k Gi k 

< P { — {Si - ^-Sk) < lOOlogz + — for 1 < z < A;/2, and —{Si - \Sk) < 
\ai k cTi ai k 

100 \og{k -i) + — for fc/2 < z < A; ) < C{1 + yf/k, 
^1 / 

where C is independent of n, k and 



(2^ 



By all the above estimates (26), (27) and (28), 

n/2 00 

P{3v e D„ such that S„ > M{n) +y,T^ = k)< C{1 + yYe'"^^^ ^ k 

k=n/A 



y\ "/t-5/2_ 

k=l 



(29) 



This can again be made as small as we wish, say at most by choosing y large enough, 
(iii). When n/2 < A; < 3n/4, we have 

P{3v G D„ such that S^ > M{n) +y,T^ = k) 
< P{3v e Bk such that S^ > M{k) + f{k) and S^. < M{i) + /(z) for 1 < z < n/2) 



< 



''El 



{5„>Af(fc)+/(fc), and s„,<A/(j)+/(i) for i<i<n/2}- 
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The above expectation is, by conditioning on = M{n) + x}, 



P{S',_^/2 > M{k) - M{n/2) + f{k) - x) ■ 



■P{Si - < f{i) -'-xioil<i< n/2) 

■psJM{n/2) + x)dx, 



(30) 



where 5* and S' are two copies of the random walks before and after time n/2, respectively, 

2 

and Ps„/2(^) is the density of Sn/2 ~ N{0, ^). 

We then estimate the three factors of the integrand separately. The first one, which is 
similar to (24), is bounded above by 



P{S',_^/, > M{k) - M{n/2) + f{k) -x)< 



_ (M(k)-M{n/i) + f(k)-xy 
f, 2(fe-n/2)<T2 



^Jk-n/2 



The second one, which is similar to (|28|), is estimated using Lemma |3 

i 



n/2^"/^ k 
The third one is simply the normal density 



Sn/2 < f{i) -jxioil<i< n/2) < C(l + 2y- xf/n. 



(31) 



Ps„,,(M(n/2)+x) 



(J (M(n/2) + 3:) 



V2 log 2 „ 



ne 



n 



(32) 



Therefore, the integral term (30) is no more than 



C{k - n/2Y^I\ 



V2 log 2 ry 



'l + 2y-xfe^ 



•J2 log 2 v/2 log 2 



dx, 



V2 log 2 



which is less than or equal to C(l + y)'^e ■'i ^{k — n/2) "^1"^ since (J2 < <Ji- 
Summing these upper bounds together, we obtain that 

3n/4 00 

^ P{3v e Dn such that > M{n) +y,T^ = k)< C(l + y)^e~'^^^ ^ A;" 

fc=n/2 fc=l 



3/2 



(33) 



This can again be made as small as we wish, say at most by choosing y large enough, 
(iv). When 3?7,/4 < A; < n, we have 

P{3v G D„ such that > M{n) +y,T^ = k) 
< P{3v e Bk such that > M{k) + f{k), and S^-. < M{i) + f{i) ioi I < t < k) 



< 



{Sv>M{k)+f{k), and s^i<M{i)+f{i), for i<j<A:}' 
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The above expectation is, by conditioning on = M{n) + x}, 



P(5U/2 > M{k) - M{n/2) + f{k) - X, 

S- < M{i) - M{n/2) + f{i) - X, for n/2 < i < k) 
■P{S^ - < /(^) - for 1 < ^ < n/2) ■ ps„„{M{n/2) + x)dx 

where S and 5" are copies of the random walks before and after time n/2, respectively. 



The second and third probabilities in the integral are already estimated in (31) and (32). 
It remains to bound the first probability. Similar to (26), it is bounded above by 



P > M{k) - M{n/2) + f{k) - X, S[ < M{i) - M{n/2) + f{i) - 

for n/2<i<k^< C{1 + 2y- xfe''^^^^""\n - k)-^/^. 



X, 



With these estimates, we obtain in this case, in the same way as in (iii), that 

n oo 

J2 Pi^v e Bn such that > M{n) + y, r„ = k) < C(l + yfe'^^^ ^ k-^'"^. 

fc=3n/4 k=l 



(34) 



This can again be made as small as we wish, say at most by choosing y large enough. 



Summing (l2§, (|29]), (|33]) and (|34j), then (|23]) and thus ([22]) follow. This concludes the 
proof of Theorem |2} □ 



4 Further Remarks 

We state several immediate generalization and open questions related to binary branching 
random walks in time inhomogeneous environments where the diffusivity of the particles 
takes more than two distinct values as a function of time and changes macroscopically. 

Results involving finitely many monotone variances can be obtained similarly to the 
results on two variances in the previous sections. Specifically, let > 2 (constant) be the 
number of inhomogeneities, {cxf > : z = 1, . . . , fc} be the set of variances and {tj > : ? = 
1, . . . , A;}, satisfying Yl\=i U = '^, denote the portions of time when af governs the diffusivity. 
Consider binary branching random walk up to time ra, where the increments over the time 
interval [Xlti ^j"-; Yl{=i ^ii^) ^(0' 1 < j < ^- When < cr| < ■ ■ ■ < are 

strictly increasing, by an argument similar to that in Section |2} the maximal displacement 
at time n, which behaves asymptotically like the maximum for independent random walks 
with effective variance Yli=i'^i^i^ 



. 2{\og2)J2Uafn-]- ^ Jll^' logn + Op(l). 
\ ^ 2 V21og2 

When af > a2 > ■ ■ ■ > <jI are strictly decreasing, by an argument similar to that in Section 
[3} the maximal displacement at time n, which behaves like the sub-maximum chosen by the 
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previous greedy strategy (see ([6])), is 



i=l ^ 



Results on other inhomogeneous environments are open and are subjects of further study. 
We only discuss some of the non rigorous intuition in the rest of this section. 

In the finitely many variances case, when {af : i = 1, . . . ,k} are not monotone in i, 
the analysis of maximal displacement could be case-by-case and a mixture of the previous 
monotone cases. The leading order term is surely a result of the optimization problem 



(12) from the large deviation. But, the second order term may depend on the fluctuation 



constraints of the path leading to the maximum, as in the monotone case. One could probably 



find hints on the fluctuation from the optimal curve solving (12). In some segments, the 
path may behave like Brownian bridge (as in the decreasing variances case), and in some 
segments, the path may behave like a random walk (as in the increasing variances case). 

In the case where the number of different variances increases as the time n increases, 
analysis seems more challenging. A special case is when the variances are decreasing, for 
example, at time < i < n the increment of the walk is A^(0, af^) with af^ = 2 — i/n. The 
heuristics (from the finitely many decreasing variances case) seem to indicate that the path 
leading to the maximum at time n cannot be left 'significantly' behind the maxima at all 
intermediate levels. This path is a 'rightmost' path. From the intuition of [11], if the allowed 
fluctuation is of order ra" [a < 1/2), then the correction term is of order n^~^", instead of 
logn in ([T]). However, the allowed fluctuation from the intermediate maxima, implicitly 
imposed by the variances, becomes complicated as the difference between the consecutive 
variances decreases to zero. A good understanding of this fluctuation may be a key to finding 
the correction term. 
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